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Abstract 

Starting from the linearized fluctuating Boussinesq equations we derive an expres- 
sion for the structure factor of fluids in stationary convection-free thermal nonequi- 
librium states, taking into account both gravity and finite-size effects. It is demon- 
strated how the combined effects of gravity and finite size causes the structure factor 
to go through a maximum value as a function of the wave number q. The appearance 
of this maximum is associated with a crossover from a q~ 4 dependence for larger q 
to a q 2 dependence for very small q. The relevance of this theoretical result for the 
interpretation of light scattering and shadowgraph experiments is elucidated. The 
relationship with studies on various aspects of the problem by other investigators is 
discussed. The paper thus provides a unified treatment for dealing with fluctuations 
in fluid layers subjected to a stationary temperature gradient regardless of the sign 
of the Rayleigh number R, provided that R is smaller than the critical value R c 
associated with the appearance of Rayleigh-Benard convection. 
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1 Introduction 



Questions concerning the nature of thermal fluctuations in fluids in thermal 
nonequilibrium states have been the subject of many studies during the past 
two decades. Specifically, in this paper we consider thermal fluctuations in a 
horizontal layer of a one-component fluid that is heated either from above or 
from below in such a way that the fluid is still in a stable convection-free state. 
A number of different and independent approaches for dealing with various 
aspects of this problem have appeared in the literature. 

Most of the earlier studies focused on the nonequilibrium fluctuations close to 
the convective Rayleigh-Benard instability because of their possible influence 
on the pattern-selection process when convection appears [1-4]. The effects of 
both linear and nonlinear terms in the hydrodynamic fluctuations close to the 
convection threshold have been evaluated. This line of research turned out to 
be very important, mainly because it lead to the development of convection 
model equations for dealing with spatiotemporal-pattern formation above the 
threshold, as reviewed by Cross and Hohenberg [5]. 

A second and independent line of research developed when investigators be- 
came interested in the nature of the fluctuations in nonequilibrium steady 
states far away from any hydrodynamic instability. The first important theo- 
retical work along this line was done by Kirkpatrick et al., who calculated the 
correct expression for the structure factor of a fluid subjected to a stationary 
temperature gradient [6] . Althought their result was originally obtained on the 
basis of the same mode coupling theory that in equilibrium leads to long-time 
tail contributions to the correlation functions for the transport coefficients, 
it turned out that the same result could also be obtained from the simpler 
scheme of Landau's fluctuating hydrodynamics [7-9]. This line of research re- 
ceived an important momentum from the experimental side when the research 
group of Sengers and coworkers verified from light-scattering experiments the 
main conclusion of the theoretical result, namely that the nonequilibrium con- 
tribution to the structure factor of a liquid is proportional to the square of the 
temperature gradient VT and inversely proportional to the fourth power of 
the wave number q of the fluctuations [10-12]. Theory and experiments have 
subsequently been extended to also deal with nonequilibrium fluctuations in 
liquid mixtures and polymer solutions, as reviewed elsewhere [13]. 

A third line of independent research has focused on the real-space behavior of 
the nonequilibrium correlations resulting from the heat-diffusion equation [14- 
17]. A study of the correlations resulting from the heat-diffusion equation is 
equivalent to a study of the nonequilibrium correlations in a liquid in the di- 
rection co-incident with the temperature gradient. This line of research has 
shown how the nonequilibrium correlation function of the fluctuating variables 
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encompasses the entire system not involving any intrinsic lenght scales. The 
real-space analysis has demonstrated the long-range nature of the nonequilib- 
rium fluctuations in the direction of the temperature gradient, which has been 
further confirmed by numerical integration of the equations [14] , by computer 
simulations [18] and by a lattice-gas automaton approach [19]. 

It is our present goal to further extend the second aforementioned research 
so as to describe nonequilibrium fluctuations far away from any convective 
instability, nonequilibrium fluctuations close to the convective instability and 
nonequilibrium fluctuations produced by the heat-diffusion equation from one 
unified point of view. The original work of Kirkpatrick et al. and of oth- 
ers [6,7,9] yielded an expression for the structure factor of a fluid in thermally 
nonequilibrium states without considering any gravity or boundary effects. 
Segre at al. extended the theory to include gravity effects, first in a one- 
component liquid [20] and then in a liquid mixture [21]. They predicted that 
gravity will cause the singular q~ A dependence of the structure factor to sat- 
urate at small values of q. This prediction was confirmed experimentally by 
Vailati and Giglio from light-scattering experiments at ultra-small scattering 
angles [22,23]. More recently, Ortiz de Zarate et al. evaluated the nonequilib- 
rium structure factor including finite-size effects but neglecting gravity [24,25]. 
They concluded that deviations from the q~ A dependence due to finite-size 
effects will be just as important as deviations due to gravity. Hence, for a 
quantitative interpretation of ultra-low-angle light scattering or shadowgraph 
experiments it is important to account for both gravity and finite-size effects. 

In this paper we extend the previous work so as to include gravity and finite- 
size effects simultaneously. Starting from the linearized fluctuating Boussinesq 
equations, we shall derive a complete expression for the structure factor of the 
nonequilibrium fluid making no other approximations than the ones incor- 
porated in the Boussinesq equations. We shall then show that the resulting 
expression is capable of describing nonequilibrium fluctuations both far away 
and close to the convective instability. Close to the convective instability we 
shall recover the linear Swift-Hohenberg correlation function as an approxi- 
mation to our more complete result. From the study of the behavior of the 
nonequilibrium correlation function in real space, we shall recover with ap- 
propriate approximations, expressions and plots previously obtained by other 
investigators for the case of the heat-diffusion equation. Thus it will be shown 
how these various conclusions can be reached from one single theoretical result. 

Attempts to arrive at a unified description of nonequilibrium fluctuations both 
far away and close to the convective instability has been made by Kirkpatrick 
and Cohen [26] and by Schmitz and Cohen [27], but the consequences of their 
theoretical results for the interpretation of actual experiments are not clear. 
In contrast, we shall be able to make a connection with available experimental 
data for the thermal fluctuations near the convective Rayleigh-Benard insta- 
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bility [28]. 



There has been an extensive debate whether fluctuating hydrodynamics can 
properly account for the thermal fluctuations close to the onset of convec- 
tion [29-33], but the most recent answer to this question appears to be posi- 
tive [28,34]. The present paper will provide additional evidence for the validity 
of fluctuating hydrodynamics to describe nonequilibrium fluctuations close to 
the onset of the classical Rayleigh-Benard instability in simple fluids. 

We shall proceed as follows. We start from the linear fluctuating Boussinesq 
equations which will be written in a form most suitable for our analysis in 
Section 2. In Section 3 we review the derivation of the well-known expression 
for the structure factor of a fluid subjected to a stationary temperature gradi- 
ent without taking into account the presence of boundaries, but including the 
effects of gravity. The relationship of the resulting expression for the effect of 
gravity with that obtained previously by Segre et al. [20] will be elucidated. In 
Section 4 we then consider the modifications to the nonequilibrium structure 
factor due to the finite height of the fluid layer. In Section 5 we present a de- 
tailed analysis of how the combined effects of gravity and finite-size will appear 
in low-angle light scattering and in shadowgraph experiments. Sections 4 and 5 
contain the most important new features of the present work. In Section 6 we 
introduce an approximation valid close to the convective instability and show 
how one can recover the linear Swift-Hohenberg model for the fluctuations 
near the convective threshold. In Section 7 we discuss the spatial behavior of 
the nonequilibrium structure factor and we shall make contact with the results 
obtained by previous investigators for the long-range nature of the nonequi- 
librium fluctuations in the direction of the temperature gradient. Our results 
are summarized in Section 8. 



2 Linearized fluctuating Boussinesq equations 

We consider a fluid layer between two horizontal plates separated by a dis- 
tance L. The fluid layer is subjected to a temperature gradient in the vertical 
direction by maintaining the plates at two different temperatures. The size of 
the system in the two horizontal X- and Y-directions is much larger than the 
size L in the vertical Z-direction. 

To determine the structure factor of the fluid we consider small fluctuations 
around the conductive solution. These small fluctuations shall be described by 
the linearized Boussinesq equations supplemented with random noise terms, 
as first considered by Zaitsev and Shliomis and by Swift and Hohenberg for 
studying the influence of thermal noise close to the convective instability [1,4]. 
Use of the Boussinesq approximation to the full hydrodynamics equations 
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implies that we neglect the sound modes and consider only density fluctuations 
caused by temperature fluctuations. The temperature gradient is applied along 
the Z-direction and is written in the form: 



VT = VT z, (1) 

where T = (T) is the average local value of the temperature T and where 
VT = dTo/dz. Hence, VT is positive when the fluid layer is heated from 
above and negative when heated from below. The gravitational force g is 
directed in the negative Z-direction: 

e = -gz, (2) 

where g is the gravitational acceleration constant. In this notation, the Rayleigh 
number may be defined as: 

_ aL 4 g ■ VTp aL^gVT 

where a is the thermal expansion coefficient, v the kinematic viscosity and 
Dt the thermal diffusivity of the fluid. 

We shall evaluate the structure factor of the fluid maintained in a convection- 
free thermal nonequilibrium state, where the average value (v) of the fluid 
velocity v will be zero. Such states correspond to both negative and positive 
values of the Rayleigh number R as long as R is smaller than the critical value 
R c . For this purpose we write the linearized fluctuating Boussinesq equations 
in the form [4,5]: 

!(V^)=,V(vS + «*(g + gW, (4a, 



= D T V 2 9-w VT + F 2 (4b) 

at 

where 6 = T — T represents the local fluctuating temperature and w is the 
fluctuating Z-component of the fluid velocity v. To eliminate the stationary 
pressure gradient from the equations we find it convenient to consider Eq. (4a) 
for V 2 u>, rather than an equation for the fluctuating fluid velocity v itself [35]. 
Finally, Fi and F 2 represent the contributions from rapidly varying short-range 
fluctuations and are related to Landau's random stress tensor ST and random 
heat flow <5Q in such a way that [36]: 

Fl = -{Vx[Vx(V-5J)]} z (5a) 
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F 2 = --IV(5Q), 



(5b) 



where p and A are the density and the thermal conductivity of the fluid, while 
the subscript z in Eq. (5a) indicates that Fi has to be identified with the 
Z-component of the vector between the curly brackets. We note that we are 
here considering only the effects of additive noise represented by F\ and F2 in 
Eqs. (4). Multiplicative noise has been also considered in the literature [37], 
but that subject lies outside the scope of the present paper. Since (9) = 
and (w) = everywhere in the fluid layer, we note that the solution of the 
fluctuating Boussinesq equations will be independent of the sign adopted for 
9orw [4,5,35,36]. 

The conditions for the validity of the Boussinesq equations, also referred to as 
Oberbeck-Boussinesq equations, have been discussed extensively in the litera- 
ture [38-40] . Specifically, in deriving the Boussinesq equations from the more 
general hydrodynamic equations, one neglects the temperature dependence of 
the various thermophysical properties, except for the density which is taken 
to vary linearly with temperature. Deviations of this so-called Boussinesq ap- 
proximation do affect the characteristic nature of pattern formation upon the 
appearance of convection [5,41-43]. However, the intensity of the nonequi- 
librium fluctuations in convection-free states below the convective instability 
appears to be much less sensitive to deviations from the Boussinesq approxi- 
mation. Segre et al. have analyzed and measured the temperature fluctuations 
in liquid toluene subjected to temperature gradients up to 220 K cur 1 [12]. 
At the higher value of the temperature gradient some of the thermophysical 
properties varied significantly over the height of the fluid layer. Nevertheless, 
the intensity of the observed nonequilibrium temperature fluctuations was 
equal to the intensity predicted when all thermophysical-property values are 
taken at their average value in the fluid layer. Similar results were obtained 
by Ahlers, measuring the critical Rayleigh number for the onset of Rayleigh- 
Benard convection and the heat transfer slightly above the threshold [42]. 

We note that in the derivation of the Boussinesq equations it is also assumed 
that the adiabatic temperature gradient (aT /cp)g, where cp is the isobaric 
specific heat capacity and T the average temperature in the fluid layer, is small 
and can be neglected compared to the magnitude of the imposed temperature 
gradient. In practice, this is a very good approximation [38,39]. In addition, 
in the Boussinesq Eq. (4b) the coefficient multiplying V 2 6> is usually identified 
with the thermal diffusivity of the fluid D T . For consistency, we have also 
expressed the prefactor of the second random noise term in Eq. (5b) in terms 
of the same diffusivity D T [39] . 
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3 Bulk structure factor of a fluid subjected to a stationary tem- 
perature gradient 



Starting from the fluctuating Boussinesq equations, we first derive an expres- 
sion for the structure factor of a fluid subjected to a stationary temperature 
gradient in the presence of gravity but in the absence of any boundary condi- 
tions. Althought the resulting expression for such a bulk structure factor will 
only be meaningful for negative Rayleigh numbers, it will enable us to develop 
a procedure most suited for subsequent incorporation of boundary effects. 

In the absence of any boundary conditions, a temporal and spatial Fourier 
transformation can be applied to Eqs. (4) to obtain a set of equations for 
the fluctuations in the vertical component of the velocity w(iv,q) and for the 
fluctuations in the temperature 6(iv, q) as a function of the frequency u and 
the wavevector q: 



V 



■q 2 {\ u + v q 2 ) a g q 2 
VT i oj + Dt 







_ fe(^q) N 


1 




\F 2 ( y iv,q) J 



(6) 



where q\\ = yq 2 + q 2 represents the magnitude of the component of the 
wavevector q in the XY-plane, i.e., the component of q perpendicular to the 
temperature gradient. We note that in this paper the parallel and perpendicu- 
lar directions are defined with respect to the horizontal plane, as done by most 
researchers in the field [8,26,27,31], but unlike the notation used by Segre et 
al. [20,21]. The Fourier transforms iq(u;,q) and F 2 (uJ,q) of the random noise 
terms can be expressed as a function of the Fourier transforms of the random 
stress tensor and the random heat flow: 

Fi (u, q) = - {g 2 [q • ST(u, q)] 2 - q z q ■ [q • 5T(u, q)] } , (7a) 



D T 

F 2 (uJ, q) = -i — q • SQ(u, q). 



(7b) 



The solution for w(cu, q) and 9(u, q) can be readily obtained by inverting the 
matrix in Eq. (6), so that: 



w(u,q) 



i iv + Dt q 2 —agq 2 



{ -VT -q 2 (i u + v q 2 )j (Fi(w,q) 



K e{u,q)J ? [w-in+(q)] [w-in_(q)] ^F 2 (u;,q)y' 



(8) 
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where we have introduced the quantities: 




(9) 



Stricly speaking, in the expression (9) for the decay rates fi±(q) the tempera- 
ture gradient VTo should be identified with the effective temperature gradient 
VTo + (aTo/cp)g, as shown by Segre et al. [20]. However, as mentioned in Sec- 
tion 2, the contribution (aT /cp)g from the adiabatic temperature gradient 
is neglected in the Boussinesq approximation. In the limit g — * 0, we have 
+ = vq 2 and = D T q 2 , which are the decay rates associated with the 
transverse- velocity fluctuations w and the temperature fluctuations 9, respec- 
tively. 

In the Boussinesq approximation, density fluctuations are caused only by tem- 
perature fluctuations, while pressure fluctuations are neglected. Consequently, 
the structure factor S(uj, q) will be related to the autocorrelation function of 
the temperature fluctuations: 



To deduce the autocorrelation function (9*(cu, q) 9(cu', q')) of the temperature 
fluctuations from Eqs. (8), we need the correlation functions for the Langevin 
noise terms Fi(u, q) and F^u, q), defined by Eqs. (5). In nonequilibrium fluc- 
tuating hydrodynamics it is assumed that the random stress tensor and the 
random heat flow retain their local-equilibrium values [8,44,45]. Thus, as in 
equilibrium, we shall consider both 5J as 5Q as "white" noise. Their autocor- 
relation functions are short ranged in time and in space, being proportional 
to delta functions. Fourier transforming the explicit expressions for the corre- 
lation functions of ST(r,t) and 5Q(r, t) as, for instance, given by Eqs. (3.12) 
in Ref. [8], and using the definition (5) of Fi(u,q) and F 2 (o;,q), we obtain: 



where k-Q represents Boltzmann's constant. Combining Eqs. (8), (10) and (11) 
we obtain for the dynamic structure factor S(u>, q) of a nonequilibrium fluid 



<0>,q) 0(o/,q')> 



1 



S(u, q) (2tt) 4 5(lu - u') 5(q - q'). (10) 



a 2 p 2 



(F*(u;,q) F l {u'^)) = 2k B T- q\ q 4 (2tt) 4 6(u - cS) 6(q-q'), 

9k T 2 \ 

(FZ( U , q) F 2 (u/, q')> = q 2 (2tt) 4 5(uj - J) 5(q - q'), 

P c p 

(F*(u,q) F 2 (u/,q')) = <^2>,q) F^u' ,q>)) = 0, 



(11) 
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in the Boussinesq approximation the following expression: 



Qf , / 2Dxq 2 (c 2 + uY) 2u(c P /T) VT 2 gf 



[cu 2 + 5^(q)] [cu 2 + ft 2 (q)] [w 2 + ^(q)] [cu 2 + fi 2 (q)] j ' 

(12) 

where Se represents the intensity of the fluctuations in thermodynamic equi- 
librium: 

Se = p 2 K T k B T?—±. (13) 

7 

Here 7 denotes the heat-capacity ratio cp/cy and k t is the isothermal com- 
pressibility. In deriving Eq. (12) we have employed the thermodynamic relation 
a 2 D T = [(7 - l)/7]A/t T /T. 

The equation for the Rayleigh spectrum of a fluid subjected to a stationary 
temperature gradient including the contribution from gravity was first derived 
by Segre et al. [20]. We recover the same result, but without the (negligible) 
contribution from the adiabatic temperature gradient. In the limit g — > 0, 
Eq. (12) reduces to the expression of the Rayleigh spectrum of a fluid subjected 
to a stationary temperature gradient, first obtained by Kirkpatrick et al. [6], 
and subsequently reproduced by other investigators [7-9] . The nonequilibrium 
structure factor, as given by Eq. (12), is anisotropic, and it depends on the 
magnitude of the horizontal component q\\ of the scattering vector directly and 
through the decay rates fi±(q). 

In this paper we shall focus our attention on the static structure factor: 

5(q) = (27r)- 1 |^5(c, q ), (14) 

which determines the total intensity of Rayleigh scattering [46]. For subsequent 
use, we find it convenient to introduce dimensionless wave numbers q = qL 
and q\\ = q\\L, where L is the finite height of the fluid layer. Upon substituting 
Eq. (9) for f2±(q) into Eq. (12) and integrating the resulting expression for 
S(cu, q), we obtain: 



S( q ) = S E ll + S° NE fzR . (15) 




In Eq. (15) S^ E represents the strength of the nonequilibrium enhancement of 
the structure factor, which is given by: 

9° - aR 1 (^-1)(cp/T)L 4 

_ cM-T ^D 2 (VT ° } ' (16) 
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where a = u/Dt is the Prandtl number. Since we have not yet considered any 
finite-size effects, Eqs. (15) and (16) do not depend explicitly on the height 
L. Equations (15) and (16) are identical to Eq. (2.35) in Ref. [20]. The term 
in Eq. (16) proportional to R is related to the adiabatic temperature gradient 
and is in practice negligibly small. Hence, the intensity of nonequilibrium 
fluctuations continues to be proportional to (VT ) 2 in the presence of gravity. 

Note that (15) can only be valid for negative Rayleigh number R, i.e., when 
the fluid layer is heated from above so that VT > 0. For any R > 0, the 
nonequilibrium contribution to the structure factor will always diverge at some 
finite value of the horizontal component q\\ of the scattering vector. The ap- 
pearance of this divergence is a direct consequence of us having performed 
the calculation without taking into account boundary conditions. As will be 
demonstrated in the next section, when boundary conditions are taken into 
account, a different expression for S(q) is obtained. This new expression will 
be valid not only for negative Rayleigh numbers, but also for a finite interval 
of positive Rayleigh numbers, up to some critical Rayleigh number R c . 



4 Modification of the nonequilibrium structure factor due to finite- 
size effects 



Since in practice the fluid layer is confined between two horizontal plates sep- 
arated by a (small) distance L, the nonequilibrium structure factor will be 
affected by the presence of boundary conditions in the Z-direction. In some 
previous publications, we have evaluated finite-size effects on the nonequi- 
librium structure factor neglecting gravity [24,25,47]. A major conclusion of 
this previous work is that finite-size effects appear at wave numbers where in 
practice the nonequilibrium structure factor will also be affected by gravity. 
Hence, for the interpretation of experimental measurements it becomes imper- 
ative to incorporate both finite-size effects and gravity effects simultaneously. 
This task is implemented in the present section. 

As in our previous publication [24], we again apply a Fourier transformation 
of the fluctuating Boussinesq equations (4) in space and in time, but restrict- 
ing the spatial Fourier transformation to the XY-plane. We thus obtain the 
following set of linear stochastic differential equations: 



/ 


' d 2 






- d 2 2 


i uj 


dz 2 ' 


- n 2 

q \\ 


— V 


dz~ 2 q \ 



VT 



i u — Dn 



a 9 q\\ 
d 2 



\ 



dz- 



-n 



W{UJ, qy, z 
6(u,q\\,z) 
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where qy is the wavevector in the XY-plane. The random noise terms F\ (u, qy , z) 
and F 2 (u;, qy, z) in Eqs. (17) are related to the partial Fourier transforms 
5T(u>, qy, z) of the random stress tensor and 5Q(cu, qy, z) of the random heat 
flux. The actual expressions are a bit complicated and can be found else- 
where [25]. 

As is often done in the literature [26,35], for the sake of simplicity, we adopt 
stress-free boundary conditions for the vertical velocity and perfectly conduct- 
ing walls for the temperature, so that: 

Q(u, qy, z) — at z — 0, L , 
w{u, qy, z) — at z — 0, L , (18) 
w(u,qu, z) = at z — 0,L. 



dz 2 



Note that these boundary conditions imply the absence of any possible fluc- 
tuations in the temperature and velocity of the fluid adjacent to the walls. As 
is well known, stress-free boundary conditions represent a fluid bounded by 
two free surfaces, which is a rather unrealistic case [35,48]. For the realistic 
case of a fluid confined between two rigid surfaces, no-slip boundary condi- 
tions in the velocity are more appropriate. For the particular case of g = 0, we 
have evaluated the nonequilibrium structure factor using both stress-free [24] 
and no-slip [25] boundary conditions. While there are numerical differences 
between the two solutions, the dependence of the nonequilibrium structure 
factor on the wave number q appears to be qualitatively the same for the two 
sets of boundary conditions. For mathematical simplicity, we evaluate here 
the nonequilibrium structure factor using stress-free boundary conditions in 
the expectation that also in the presence of gravity the more realistic no-slip 
boundary conditions will not modify the qualitative nature of the wave- number 
dependence of the structure factor. 

To search for a solution of Eq. (17), we represent w(u,q\\,z) and 9(u,q\\,z) 
as a series expansion in a complete set of eigenfunctions of the differential 
operator in Eq. (17), satisfying the boundary conditions (18) [24]. Because of 
the simplicity of the boundary conditions considered here, an appropriate set 
of eigenfunctions is the Fourier sine basis in the [0, L\ interval [35]. We thus 
assume: 



^w(u, qy, z 
K d(u, q\\,z) J aCi \B n (u, qy) J 



{Nnz\ . . 

sin — . (19) 



To deduce the coefficients A N (u,ci\\) and B N (u,ci\\) from Eq. (17), we also 
have to represent the random noise terms Fi(u, qy, z) and F 2 (u;,qy,z) as a 
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Fourier sine series: 



F 1 (u,q h z) \ _ ~ / F ljN (u, qy ) j ^ /Nirz" 
F 2 (uj,q h z) N=i\F 2>N (u,q^' ^ 1 



(20) 



where we have introduced the set of random functions Fi jN (u, qy) and F 2>N (ui, qy), 
which are the projections of the random noise terms over the eigenfunction 
basis. They are given by: 



Fi )JV (a;,qy 
\F 2 ,n(u, qy 



2 r L 



Fi(w,q||,z) \ , /JVtiV. , 
] sin I — - — I ctz. 

F 2 (a;,q||,z) 



(21) 



Representing the random noise terms by Eq. (20), one readily deduces from 
Eq. (17) expressions for the coefficients of the Fourier series A N (u,q\\) and 
5jv(w,qy). To calculate the structure factor in the Boussinesq approximation, 
we only need the coefficients B^(u, qy) for the temperature fluctuations. Sub- 
stituting Eqs. (19) and (20) into Eq. (17) we obtain: 



-VT F ltN (u>, qy) - (gg + gj) [i u + v (ffi + gf)] i^co, qy) 



(22) 



where, similarly to Eq. (9), we have introduced decay rates Q±(N,q\\), which 
are now given by: 



[y - £>t) j 



(23) 



Since the differential operator in Eq. (17) depends only on the modulus of qy, 
the problem has cylindrical symmetry and the solution will only depend on 
the magnitude q\\ of the two-dimensional wavevector qy. 

In analogy to Eq. (10), the relationship between the dynamic structure factor 
S(u>, q\\, z, z') and the autocorrelation function of the temperature fluctuations 
is now given by: 

(9*(u,q h z) 9(u',^,z')) = ^ S(u>,q h z,2>) 5(u - u f ) 5(q l} - qj,) (24) 



To obtain the structure factor S(u, gy, z, z') we need the correlation functions 
between the different projections of the random noise terms. These have been 
calculated in a previous publication with the result [24]: 
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,1/2 2 / 2 TVV 



'NM 



(Fl N (u, q,|) • F liM (u/, qj,)) = gg + 

x(2tt) 3 <J(q||-qJ|), 

(Fl N (u, q,|) • F 2 , M (u/, qj,)) = (F 2 >(o;, q,|) • F 1>M (u/, qj,)) = 0, (25) 

/™ / x r, / / ,u 2/c B T 2 A 2 / 2 iV 2 7r 2 \ p 
(^2,i\r(^q||)-- F 2,Af(^,q||))= ^ l g ll + _ L^J ^ M 

x (2tt) 3 % -u/) <J(q|| - qj,). 

As in the case of the absence of any boundary conditions, we continue to 
assume that the correlation functions between the different components of 
the random current tensor and the random heat flux retain their equilibrium 
values. This assumption remains valid as long as L is a macroscopic distance, 
much larger than the molecular distances in the fluid. 

In this paper we are interested in the static structure factor S(q\\,z,z') = 
(27r) _1 / duS(uj, q», z, z'). With the information presented above, this quantity 
can be readily calculated. In analogy to Eq. (15) for S(q) in the absence of 
any boundary conditions, the modified nonequilibrium structure factor can be 
written in the form: 

S(q\\,z,z) = S E [d(z - z') + S^ E S NE (q\\,z,z')] , (26) 

where Se and S^ E are again given by Eqs. (13) and (16), respectively, while 
5ne(<Z||, z, z') incorporates the finite-size effects and is given by: 

{Nixz\ (Nixz'\ 



The first term in Eq. (26) is again the static structure factor of a fluid in ther- 
modynamic equilibrium, which is not affected by any finite-size effects [46]. 
The second term in Eq. (26) represents the nonequilibrium enhancement of the 
structure factor. This nonequilibrium enhancement is proportional to (VTo) 2 
through the expression (16) for S^ E , it depends on the gravitational accelera- 
tion constant g through the appearance of the Rayleigh number in Eqs. (16) 
and (26), and it depends on the finite height L of the fluid layer explicitly 
in Eq. (26) and also through qu = q»L. In the limit R — > we recover the 
expression for S(q\\,z,z') presented in a previous publication [24]. 

It is important to note that Eq. (27) is valid for both negative and positive 
Rayleigh numbers, provided that 

R<R C = ^. (28) 
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For R> R c there always exist values of q\\ for which right-hand side of Eq. (27) 
diverges. Of course, the value R c = 277r 4 /4 equals the well-known value ob- 
tained form a linear stability analysis of the Boussinesq equations with no- 
slip boundary conditions [35]. As expected, for R > R c the present calcula- 
tion breaks down because the fluctuations around the conductive state are no 
longer small. 



Equation (27) for the normalized nonequilibrium enhancement Sne(<Z||, z, z') 
represents our principal result for the combined effects of gravity and finite- 
size on the nonequilibrium structure factor of the fluid. The remaining part of 
this paper is concerned with an analysis of some of the physical consequences 
that follow from Eq. (27). 

The double Fourier series contained in Eq. (27), when it is convergent, can be 
easily summed. To do so, we first find the roots of the denominator for iV 2 . 
We then split the inverse of the product of three binomials as a sum of three 
terms, and we finally obtain for Sne(<7||, z, z'): 



sin 



\,Z,Z 



i?*T E E 



Nirz 
L 



sin 



^ Nirz r 



3L 



N=l i=0 



A? 



N 2 7l 2 



4/3 



- 1 



(29) 



where A; for % — 0, 1, 2 are the three complex cubic roots of R: 



A = IRI 1 / 3 



expi | 
expi^ 



yexpi i£±*2L 



(30) 



In Eq. (30) represents the phase of R when considered as a complex number: 
ip = for R > and ip = ir for R < 0. This procedure is only valid for R ^ 0. 
However, the finite-size effects on the structure factor for the particular case 
R = have already been considered in a previous publication [24], where the 
sum of the series in Eq. (27) for this particular case has been performed by 
a different method. Hence, we do not need to re-evaluate Sne for the case 
R = 0. To implement the summation in Eq. (29) we consider the relation: 



E 



N= 



[Ntxz\ 
sin I — \ sin 



— 7T 



C0S r#(L 




(31) 



co S r# ( L 



\z + z'\) 



Ayfji sin (iry/Ji) 

This equation can be obtained from formula 1.445 in Ref. [49] and it is valid 
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for z, z' G [0, L\ and for any complex number /x, provided that in the right- 
hand side of (31) the square root of /i with positive real part is chosen. Now, 
the sum of the double Fourier series contained in Eq. (27) can be performed, 
but the result is a long expression, not particularly informative, which can be 
easily obtained by the reader. Explicit expressions to be obtained upon further 
integration of Sne(<Z, z, z') will be presented in the following section. 



5 Consequences for light scattering and shadowgraph experiments 

5.1 Nonequilibrium structure factor probed in experiments 

The nonequilibrium fluctuations can be obserbed in small-angle light-scattering 
experiments. For a discussion of such light-scattering experiments we con- 
sider an experimental arrangement like the ones employed at the University 
of Maryland [9-12] and by Vailati and Giglio at the University of Milan [22,23]. 
A schematic representation of such a light-scattering experiment is shown in 
Fig. 1. The scattering medium is a thin horizontal fluid layer bounded by two 
parallel plates whose temperatures can be controlled independently so as to 
establish a temperature gradient across the fluid layer. The temperature gra- 
dient can be parallel or antiparallel to the direction of gravity. The horizontal 
plates are furnished with windows allowing laser light to propagate through 
the fluid in the direction parallel to the gravity and to the temperature gra- 
dient. Light scattered over an angle <fi arises from fluctuations with a wave 
number such that [46]: 

g = 2g o sin(0/2), (32) 

where go is the wave number of the incident light inside the scattering medium. 
To observe any nonequilibrium fluctuations one needs to observe the scattered 
light at small wave numbers and, hence, at small scattering angles. 

From electromagnetic theory [46] it follows that the scattering intensity S'(q) 
is obtained from an integration of the structure factor over the scattering 
volume, so that [24,27]: 

L L 

S(q h q±) = lf J e-^-*') S(q h z,z>) dz dz'. (33) 



In Eq. (33) we have assumed that the scattering volume extends over the full 
height of the fluid layer as is the case in small-angle light scattering from a 
thin fluid layer [12]. In this situation scattered light received in the collecting 
pinhole of the detector indeed arises from all points illuminated by the laser 
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Fig. 1. Schematic representation of a nonequilibrim light-scattering experiment, qi 
is the wave vector of the incident light and q s is the wave vector of the scattered 
light. The magnitude q = |q; — q s | of the scattering wave vector is related to the 
scattering angle <f> by q = 2qo sin(0/2), where go is the magnitude of the wave vector 
qi of the incident light inside the fluid. 

beam inside the fluid layer. From Eq. (32) and the geometrical arrangement 
shown in Fig. 1, we note that q\\ and q± in an actual light-scattering experiment 
are not independent variables, because they are related to the scattering angle, 
0, by: 



Equation (34) shows that for actual small-angle light-scattering experiments 
one may use the approximation q\\ ~ q, q± ~ 0. 

An alternative promising experimental technique for measuring the intensity, 
though not the temporal spectrum, of nonequilibrium fluctuations is quan- 
titative shadowgraph analysis [28,50-53]. The experimental arrangement is 
similar to the one depicted in Fig. 1 but, instead of a laser beam, an extended 
uniform monochromatic light source is employed to illuminate the sample. 
Then many shadowgraph images of a plane perpendicular to the temperature 
gradient are obtained with a CCD detector, which measures the spatial distri- 
bution of intensity /(x), where x is a two-dimensional position vector in the 
imaging plane. For each image a shadowgraph signal is defined by: 



where /o( x ) is the blank intensity distribution, when there are no fluctuations 
in the index of refraction of the sample. In practice /o( x ) is calculated by av- 
eraging over many shadowgraph images, so that fluctuations cancel out and 



q\\ = q cos (0/2) = 2g sin (0/2) cos (0/2), 
q ± = q sin (0/2) = 2q sin 2 (0/2). 



(34) 




(35) 
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the resulting lo( x ) contains only contributions from nonuniform illumination 
of the sample. From physical and geometrical optics one can demonstrate that 
the modulus squared of the two-dimensional Fourier transform of the shad- 
owgraph signal, |i(q)| 2 , after taking an azimuthal average, can be expressed 
as [50,51,53]: 



\*(<l)\ 2 (q)=T(q) S(q {l =q,q ± = 0), (36) 

where the overline indicates the azimuthal average, in which case the result 
depends only on the modulus q\\ of the two-dimensional Fourier vector q. The 
factor T(q) represents an optical transfer function, which can be derived from 
the optical arrangement used to produce the shadowgraph pictures; it includes 
contributions from the response of the CCD detector and the density depen- 
dence of the refractive index [50,51,53]. Therefore, applying a two-dimensional 
Fourier transform to the shadowgraph images, one can deduce the structure 
factor of the fluid as a function of the wave number q at q±_ = 0. 

There exist an equivalence between small-angle light scattering and shadow- 
graphy, in the sense that both methods give us S(q\\ — q, q± — 0). For light 
scattering, q is the scattering wavevector as given by (32), whereas for shad- 
owgraph techniques q is the modulus of the two-dimensional Fourier vector in 
the imaging plane. As mentioned by Bodenschatz et al., what it is measured 
in the experiments is a kind of vertical average of the fluctuations [51]. 

To obtain an explicit expression for the structure factor S(q\\,q±), we substi- 
tute Eq. (26) with Sne(9||, z , z> ) given by Eq. (27) into Eq. (33), implement 
the summation by using the relation (31) and perform the double integra- 
tion in Eq. (33). Introducing again dimensionless wave numbers q\\ = q\\L and 
q± = q±L, we thus obtain: 

S(q h q±) = S E {l + S° E S NE (q h q ± )} (37) 

with 



q — nq« 



n 2 



cos (q ± ) - 


- cos 


( 9 n 


/a,- 




sin ^ 




jx 3 - 







3*- - A,,f 3 

In Eq. (37), Se and S^e are a g & i n given by Eqs. (13) and (16), respectively, 
while Aj in Eq. (38) again represent the three complex cubic roots of R de- 
fined by Eq. (30). In thermodynamic equilibrium VT = and R = 0, so 
that S(q\\,q±) = S E and we recover from Eq. (37) the well-known isotropic 
equilibrium structure factor given by Eq. (13). 
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If we retain only the first term on the right-hand side of Eq. (38) and substitute 
it in Eq. (37) we recover the expression (15) for the bulk structure factor 
of the fluid. As noted in Section 3, this expression for the bulk structure 
factor breaks down for any positive value of the Rayleigh number because 
of a divergence at q 6 = q 2 RL~ 4 . The second term in Eq. (38) for £ne(?||,?±), 
containing a sum over the three complex cubic roots of R, represents the finite- 
size effects on the nonequilibrium structure factor and depends on the height 
L of the fluid layer. Although perhaps not immediately transparent, it can be 
verified that this second term when multiplied by S^ E vanishes in the limit 
L — > oo. Moreover, for any finite value of L, the divergence of the first term 
in Eq. (38) at q 6 = q 2 RL~ 4 for positive R is now exactly compensated by a 
similar divergence in the second term of Eq. (38) for the same wave number, 
since A = R 1 ^ 3 for positive R. Thus Sne(<7||, 3±), as given by Eq. (38), is 
continuous at q 6 = q 2 RL~ 4 for positive R. However, due to the presence of a 
sine term in the denominator, new divergences appear when: 



qWyl = Nir (39) 

for A" = 0, 1, ... . For A^ = this new divergence is again compensated due 

to the presence of the term i? 1 / 3 — q^ 3 in the numerator. But for N — 1, 

there is no compensation and the contribution (38) to the nonequilibrium 
enhancement of the structure factor diverges. Since condition (39) for A^ = 1 
is satisfied when R = (q 2 +rc 2 ) 3 /q 2 , this divergence will appear for any Rayleigh 
number larger than the critical Rayleigh number R c = 277r 4 /4. For R = R c , 
condition (39) is satisfied at a critical wave number qn c = ix j 'y/2. We thus 
recover from an analysis of the divergences in the nonequilibrium structure 
factor the well-known results of the linear stability theory for a fluid layer 
bounded by two free surfaces [35]. A more detailed analysis of the behavior 
of our solution for the nonequilibrium structure factor for positive R near R c 
will be presented in Section 6. 

As mentioned earlier, for the interpretation of low-angle scattering or shad- 
owgraph experiments one may consider the approximation q\\ ps q, q± ps 0. In 
the limit q± — > 0, Eq. (38) for Sne(9||, <?_l) can be slightly simplified; in this 
approximation the nonequilibrium enhancement Sne(<?) will depend only on q, 
which is the experimentally relevant quantity. In the remainder of this paper 
the small-angle approximation will always be assumed. We note that all the 
divergent features discussed above remain the same in this approximation. 

It is interesting to look at the asymptotic behavior of the nonequilibrium 
enhancement for q — > and for q — > oo. From Eq. (38) it can be readily shown 
that: 

^NE(g) ^ <f , (40) 

w 20160 H ' V ; 
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and 



SMQ) ^ j,- (41) 

We note from Eqs. (40) and (41) that the limiting behavior of Sne(<?) for q — > 
and q — > oo does not depend on the value of the Rayleigh number R. 

5.2 Interpretation of experimental results 

In accordance with Eq. (37), the product S^ E Sne(<?) represents the nonequi- 
librium enhancement of the structure factor with S^ E given by Eq. (16) and 
Sne(?) given by Eq. (38) with q± — > 0. From Eq. (16) we see that the pref- 
actor S^e depends on various thermophysical properties of the fluid and on 
the magnitude of the applied temperature gradient, but S^ E is independent 
of the wave number q. Hence, to evaluate the wave-number dependence of the 
nonequilibrium enhancement of the structure factor, we only need to focus our 
attention on Sne(<?)- 

Assuming q± = 0, we have evaluated 5ne(<?) from Eq. (38). In Fig. 2 we show 
on a double-logarithmic scale Sne^) as a function of q for three different val- 
ues of the Rayleigh number, namely for a large negative value, R = —5000, 
for a value R ~ which corresponds to the case investigated in a previous 
publication [24] and for a value R = 600 which is close to the critical value 
R c ~ 656. At larger values of q, Sne(<?) varies as g~ 4 , independent of the 
Rayleigh number, in agreement with Eq. (41). Upon decrease of q, Sne(<?) 
goes through a maximum and for very small values of q, S]<se{q) decreases as 
q 2 , in agreement with Eq. (40). For positive R, Sne(qD develops a prominent 
peak close to q c = n/y/2, which diverges as R — > R c , as further discussed 
in Section 6.1. Hence, a major effect of the additive noise terms in the fluc- 
tuating Boussinesq equations is the appearance of (fluctuating) patterns in 
the fluid, even below the convective instability as discussed by some other 
investigators [54,55]. 

Sengers and coworkers have measured the nonequilibrium fluctuations in liquid 
toluene [10-12] and in liquid n-hexane [56]. These experiments correspond to 
Rayleigh numbers from —25, 000 to —300, 000 at dimensionless wave numbers 
ranging from q = 640 down to q — 345. The experiments have provided an 
accurate confirmation of the q~ 4 dependence of the intensity of nonequilibrium 
fluctuations in this range of wave numbers. 

Giglio and coworkers have measured the intensity of nonequilibrium fluctua- 
tions for negative Rayleigh numbers down to wave numbers q of order unity, 
both with ultra-low angle light-scattering experiments [22,23] and from quanti- 
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Fig. 2. Normalized nonequilibrium enhancement Sne(q) of the structure factor as 
given by Eq. (38), as a function of the dimensionless wave number q for three 
different values of the Rayleigh number R. The solid curve is for R = 600, which is 
near the convective instability, the dashed curve is for R « and the dotted curve 
for R = -5000. 

tative analysis of shadowgraph images [52,53]. They actually measured the in- 
tensity of nonequilibrium concentration fluctuations in a liquid mixture. How- 
ever, due to the similar structure of the underlying hydrodynamic equations, 
the q dependence of the contribution of nonequilibrium concentration fluctua- 
tions to the structure factor in a liquid mixture is expected to be similar to the 
q dependence of the contribution of nonequilibrium temperature fluctuations 
to the structure factor of a one-component fluid [13,47]. Giglio and coworkers 
have not only confirmed the q~ 4 dependence of the nonequilibrium structure 
factor, but they have also observed the crossover to a region of q close to unity 
where the nonequilibrium structure factor is independent of q, in agreement 
with the flat region indicated in Fig. 2 for Sne(<?) at large negative Rayleigh 
numbers. 

Wu et al. have measured the nonequilibrium structure factor in a layer of fluid 
carbon dioxide, at a pressure of about 3 MPa, near the convective instabil- 
ity. For positive Rayleigh numbers, 5ne(<?) strongly depends on the parameter 
e = (R — R c ) / R c which measures the distance from the Rayleigh- Benard in- 
stability. In Fig. 3 we show Sne(<?) as a function of q close to q c = 
as calculated from Eq. (38) for e = -0.03, e = -0.01 and e = -0.005. We 
did not attempt to make a quantitative comparison with the measurements 
of Wu et al., because the slip- free boundary conditions adopted in this pa- 
per do not correspond to the experimental boundary conditions. However, a 
comparison of Fig. 3 with Fig. 3 in Ref. [28] shows that our solution yields 
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Fig. 3. Normalized nonequilibrium enhancement Sne(q) of the structure factor as 
given by Eq. (38), as a function of the dimensionless wave number q for three 
different values of the distance e to the convective instability. The solid curve is 
for e = — 5 x 10~ 3 , the dotted curve is for e = — 10~ 2 and the dashed curve for 
e = -3 x 10~ 2 . 

a good qualitative representation of the q dependence of the nonequilibrium 
structure factor as measured by Wu et al. We conclude that our solution of the 
linearized fluctuating Boussinesq equations for the nonequilibrium structure 
factor is consistent with the characteristic features of the nonequilibrium struc- 
ture factor observed in experiments for both negative and positive Rayleigh 
numbers. 



6 Nonequilibrium fluctuations close to the convective instability 

The nature of thermal noise near the convective instability has been the sub- 
ject of studies by many investigators [1,4,26,27,31,36]. Hence, it is of interest 
to make a comparison of those results with our solution for the intensity of 
temperature fluctuations for thermal nonequilibrium states. It turns out that 
for an analysis of Sne(?) near the convective instability it is more convenient 
not to perform the summation of the series in Eq. (27). Instead, substituting 
Eq. (27) into Eq. (33) and integrating the terms in the series individually, we 
can write Sne(<?) as: 



S NE (g) = 4 ]T 



oc 



q 2 1 - cos(iV7r) 



(42) 



(<f + iV% 2 ) 3 - R q 2 N 2 n 2 



which is equivalent to the expression (38). 
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6. 1 Divergence of nonequilibrium fluctuations at the convective instability 



Zaitsev and Shliomis [1] were the first to compute thermal fluctuations in a 
fluid layer subjected to a stationary temperature gradient near the convective 
instability. Using linear perturbation theory they found that the structure 
factor diverges as (R c — The same divergence follows from our solu- 

tion for Sne(<7)- To reproduce the divergence we first calculate the value q max 
corresponding to the maximum of Sne(<?) by expanding g max in powers of 
e — (R — R c )/R c around rr/^/2. Differentiating the expression (42) for Sne(<?) 
and demanding the derivative to vanish, we find that the position of the max- 
imum corresponds to: 



vr 8lV2n " (iV 2 -l)(l + 2iV 2 ) 2 (l-cosiV 7 r) 3 

W-^ + e g 2. N2 [(1 + 2iy2)3 _ 27] 2 



^ ^ + 6.324 x 1(T 4 e 2 . (43) 
V2 



Substituting Eq. (43) into Eq. (42) we conclude that the structure factor, 
which is proportional to the intensity of the scattered light, diverges when the 
convective instability is approached such that: 

~ . , -547T 2 16^ l-cosiV7r ^. , , A A . 

^ne(w) = — + ~ 6 E iv 2[(1 + 2i v 2 )3_ 27] + °^ ( 44 ) 

We thus recover the linear divere; ence of 5*ne &s a function of (R — R c ) 1 
obtained by Zaitsev and Shilomis [1] and confirmed by Swift and Hohen- 
berg [4,36]. Linear fluctuation theory amounts to a mean-field theory of fluc- 
tuations. Extremely close to the instability nonlinear effects will cause a smear- 
ing out of the sharp mean- field transition [2,4,36], but this effect will only be 
noticeable for very small values of |e| ~ 2.9 x 10~ 5 [40]. Hence, observation 
of this linear divergence is possible in experiments [28]. Deviations from lin- 
ear fluctuation theory have been observed by Scherer et al. in the case of 
electro convection [57]. 



6.2 The most unstable mode and the Swift- Hohenberg approximation 



The approximation scheme used by Zaitsev and Shliomis and by Swift and 
Hohenberg is equivalent to retaining only the term iV = 1 in the series expan- 
sion (42) for Sne(<?)- Note that, for N — 1, when R is close to R c and q\\ is 
close to q c , the denominator in the corresponding term in Eq. (27) approaches 
zero. Therefore, in this situation the term for iV = 1 is much larger than the 
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terms for any other value of N. Consequently, when R ~ R c and q ~ q c , trun- 
cating the series (42) at N — 1 yields a very good approximation. We thus 
deduce from our solution for the structure factor as measured in low-angle 
light scattering or shadowgraph experiments: 

%) = g + lr y ■ < 45 > 
— R 

where the equilibrium contribution to S(q) has been neglected ( cf., Eq. (37)). 
This is usually called the most-unstable-mode approximation in the litera- 
ture [1,4]. The denominator in Eq. (45) is zero at R — R c — 27n 4 /A and 
Q\\ = Q\\c — n /V%- Expanding the denominator in powers of gjj around gjj c , one 
obtains: 

s(q) = I* 3 <? (46) 



where, following the notation of Hohenberg and Swift [36] , we have introduced 
the parameter £q = 2/ (^/Sir 2 ). Equation (46) is proportional to the correlation 
function of the Swift-Hohenberg (SH) model below the convective instability, 
when nonlinear terms are negligible [4,36]. The results of SH refer to the 
autocorrelation of a field ip(r\\,t), which is a generalized two-dimensional or- 
der parameter introduced by them to describe the convective instability as a 
nonequilibrium phase transition. If we compare the term N = 1 in Eq. (27) 
with Eq. (2.10) in Ref. [36], we identify our Sne(<?) as proportional to the 
equal-time autocorrelation function of the field ip. Therefore the SH model 
can be obtained as an approximation close to the convective instability from 
our exact results for S(q\\ — q, q± — 0), as given by Eqs. (37) and (38). As dis- 
cussed by SH, Eq. (46) is valid for both "stress-free" and "no-slip" boundary 
conditions, but the numerical values of the constants R c , q\\ c and £q do depend 
on the boundary conditions (see Table I in Ref. [36]). An interesting feature 
of the SH model is that is a "Lagrangian" model, because the autocorrelation 
function (46) can be obtained from a simple probability functional for the field 
^(ry), as is discussed in detail in the relevant literature [4,36]. 

To compare our more complete solution of the fluctuating Boussinesq equa- 
tions for Sne(<?) with the various approximations proposed for the structure 
factor close to the convective instability, we have plotted in Fig. 4 Sne(<?) 
as a function of q. The solid line represents Sne(<?) as given by our exact 
solution (38), the dotted curve (hardly visible) the most-unstable-mode ap- 
proximation (45) and the dashed curve the SH-model approximation (46). All 
three curves correspond to a fluid layer subjected to a temperature gradient 
such that e = —0.01. From Fig. 4 it is evident that Eq. (45) yields a very 
good approximation although at large q deviations appear. The reason is that 
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Fig. 4. Normalized nonequilibrium enhancement <Sne(<7) of the structure factor as 
a function of the dimensionless wave number q near the convective instability, 
e = — 10~ 2 . The solid line represents the exact expression as given by Eq. (38). 
The dotted curve (indistingishable from the former, except asymptotically for large 
q) represents the most-unstable-mode approximation given by Eq. (45). The dashed 
curve represents the Swift-Hohenberg approximation given by Eq. (46). 

5ne(<?) for large q should vary as g~ 4 , while Eq. (45) yields (8/7r 2 )g~ 4 for large 

Many authors have represented the q dependence of the structure factor near 
the instability in terms of a Lorentzian profile centered at q c with a width 
proportional to e [26,27,31,36], and Wu et al. have also analyzed their ex- 
perimental data in terms of a Lorentzian profile [28]. However, our actual 
solution (38) does not yield a Lorentzian profile, although it can be obtained 
by expanding the denominator in Eq. (45) not in powers of q 2 around q 2 , but 
in powers of q around q c . The problem is that a Lorentzian does not lead to 
the proper asymptotic behavior for either small q or large q. It even leads to 
an apparent divergence when one tries to calculate the power of the fluctua- 
tions by integrating S(q) over all two-dimensional wave vectors q [28,36]. We 
conclude that it would be much better to represent the experimental data in 
terms of an equation like Eq. (45), which is not a Lorentzian. 

6. 3 Power of thermal fluctuations 

Another feature of nonequilibrium fluctuations close to the convective insta- 
bility that has been studied theoretically and experimentally is the so-called 
power of the thermal fluctuations, (ST 2 ). The mean-square amplitude (ST 2 ) of 
the temperature fluctuations is related to the integral of the structure factor 
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measured in the experiments: 



^-WlI s( " ] W?' (47) 

which corresponds to a vertical average of the real-space temperature auto- 
correlation function. Due to the symmetry of the problem this quantity does 
not depend on the point xy in the horizontal plane at which is evaluated. In 
Eq. (47) the prefactor l/(a 2 p 2 ) accounts for the conversion from density fluc- 
tuations to temperature fluctuations (cf., Eqs. (10) and (26)) and the factor 
XjL results from the vertical average (compare Eq. (33) with Eq. (A14) in 
Ref. [36]). From Eq. (42) we note that in terms of our solution of the fluctu- 
ating Boussinesq equations (neglecting the equilibrium contribution): 

2 S E S° NE " q 2 2ixq dq 1 - cos(iV7r) 

[ ' a 2 p 2 Lir 2 ^ Jo (q 2 + NH 2 f - R q 2 N 2 ir 2 ' 1 ' 

The integral in (48) can be performed analytically, but the result is long and 
not particularly interesting. Asymptotically close to the convective instability 
we obtain: 

(ST 2 ) ~ SeS ^ e ^ (49) 



Thus the mean-square amplitude (ST 2 ) calculated from our solution will di- 
verge as l/y/— e, which is the kind of divergence observed experimentally by 
Ahlers and coworkers [28,51]. In Eq. (49), S^ E has to be evaluated at the 
critical temperature gradient. It is interesting to note that the SH approxima- 
tion gives for (ST 2 ) the same asymptotic behavior, Eq. (49), than the exact 
solution. 



7 Correlations in real space 



In this section we study the real-space behavior of the autocorrelation func- 
tion of the nonequilibrium density fluctuations. We also shall make contact 
with the third approach dealing with nonequilibrium fluctuations that was 
mentioned in the Introduction. To calculate the autocorrelation function in 
real space we apply an inverse Fourier transform in the horizontal plane to 
the structure factor as given by Eq. (26). The equilibrium contribution gives 
a triple delta function. The nonequilibrium contribution gives a nontrivial 
contribution, which we denote as G-?m{r\\,z,z') and which is equal to: 

Se S° f°° 

G NE (r h z,z') = ^ E q\\J (q\\r\\) S NE (q h z : z') dq h (50) 
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Fig. 5. Normalized dimensionless nonequilibrium correlation function Gne^Nj z, z') 
in the Z-direction (ry = 0) as a function of z/L for z' = L/2. The solid curve is for 
R = 500, the dotted lines are for R = and the dashed curve is for R = —2000. 

where Jo(x) is the Bessel function of the first kind and order zero. Note that 
due to the cylindrical symmetry of the problem the correlation in real space 
depends only on the cylindrical radial coordinate (ry). It is difficult to study 
the integral (50) in the general case, but there are a couple of particular cases 
for which integral (50) can be evaluated explicitly or where simple approxima- 
tions can be made so that the real-space behavior of the correlation function 
can be studied in detail. 



7.1 Correlations in the vertical direction 



The first interesting case for which the long-range part of the real-space cor- 
relation function G^{r\^z,z') can be studied is for ry = 0, i.e., along the 
vertical Z-direction. In this particular case, the correlation function Gne will 
only depend on the variables z and z' and on the Rayleigh number R. Close to 
the convective instability, the behavior of the real-space correlation function 
will be dominated by the most unstable mode. From Eq. (27) we note that 
Gne w iH be proportional to the product sin^nz/L) sm(jxz' / L). Therefore, the 
concavity of Gne(0, z, z') as a function of z for fixed value of z' is negative. 
For the interesting case R = 0, Gne(0, z, z') was evaluated exactly in a previ- 
ous publication [24]. It turns out that the concavity of Gne(0, z, z') changes 
at R = 0, being negative for positive R and positive for negative R (except 
for a small region near z — z'). In Fig. 5 we have plotted the dimensionless 
nonequilibrium correlation function, G NE = G ne L 3 /(S'e5'ne)) as a function of 
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z, for z' = L/2, r\\ = and for three values of the Rayleigh number R. The 
dotted lines correspond to R = 0, the solid curve corresponds to a positive 
R and the dashed curve to a negative R. The change of the concavity of the 
function at R = is evident. This feature can be investigated analytically 
since, for small values of the Rayleigh number, Gne(0, z, z') can be readily 
expanded in powers of R yielding: 



Gne(0, z, z') 
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(51) 



where Bq(x) is the Bernoulli polynomial of the 6th degree (actually, it is possi- 
ble to obtain a systematic expansion of Gne(0, z, z') in terms of the Bernoulli 
polynomials, but the resulting series is not very informative). For R = 0, we 
recover from Eq. (51) the expression obtained for the spatial dependence of 
the correlation function in a previous publication (see Eq. (25) and Fig 2 in 
Ref. [24]). The change of concavity of Gne(0, z, z') at R = is also evident 
from Eq. (51). We note that for R = the function Gne(0, z, z') has a kink at 
z = z', indicating the presence of a discontinuity in the derivative at this point. 
For any other value of R < R c the function and its derivative are continuous. 

From Fig 5 we also note that the nonequilibrium real-space equal-time cor- 
relation function is always positive, has a maximum at z — z' and decreases 
monotonically from this maximum reaching zero at both ends of the [0, L] 
interval. Thus, we conclude that the real-space correlation function in the 
vertical direction is nonlocal, long ranged and does not involve any intrinsic 
length scale, i.e., the correlation encompasses the entire system only to be cut 
off by the finite size of the system itself. Hence, the presence of the tempera- 
ture gradient induces long-range correlations in the direction coincident with 
the temperature gradient, as has been elucidated by previous investigators, 
mainly on the basis of the nonequilibrium heat-diffusion equation [14,58,59]. 
Considering only the heat-diffusion equation means that one neglects any cou- 
pling between temperature and velocity fluctuations. This approach neglects 
also any gravity effects and is restricted to fluctuations with q\\ = 0. The 
existence of long-range correlations in the direction of the temperature gra- 
dient has been verified by numerical integration [17], by computer simula- 
tion [18], by lattice-gas-automaton dynamics [19] and by a master-equation 
approach [60]. Our result for the correlation function along the direction of 
the temperature gradient is in agreement with that obtained by Malek Man- 
sour and coworkers [14,17], as is also evident form a comparison of our Fig. 5 
for R = with Fig. 1 in Ref. [17]. Some of these authors have also con- 
sidered a coupling between temperature and velocity fluctuations [61]. But 
to simplify the hydrodynamic equations they made the assumption that the 
thermal expansion coefficient of the fluid vanishes. This is equivalent to assum- 
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r/L 

Fig. 6. Normalized dimensionless nonequilibrium correlation function Gne(^||, z, z') 
in the central plane of the liquid layer {z = z' = L/2), near the convective instability, 
as a function of r\\/L, divided by its value at r\\ = 0. The solid curve is for e = — 1CT 2 ; 
the dashed curve is the asymptotic approximation, Eq. (53), for the same value of 
e. The dotted curve is for e = —0.7. 



ing that density fluctuations are caused only by pressure fluctuations. Hence, 
their calculations are concerned with Brillouin scattering, while here we are 
considering Rayleigh scattering, for which accurate experimental information 
on nonequilibrium fluctuations has been obtained. More recently, Bena et al. 
have studied the velocity autocorrelation function adopting an hydrodynamic 
approximation known as Kolmogorov flow [62]; again, due to the nature of 
the hydrodynamic simplifications performed, the results cannot be compared 
with those obtained in the present paper. 



7.2 Correlations in the horizontal plane 



A second particular case for which the real-space behavior of the correlations 
in a liquid subjected to a stationary temperature gradient can be easily studied 
is in a horizontal plane, that is, in the plane perpendicular to the temperature 
gradient and far from the boundaries (z = z' = L/2). We shall analyze here 
the spatial dependence of the correlations in the horizontal direction near the 
convective instability using the most-unstable-mode approximation, Eq. (45), 
which, as was shown in Section 6, yields an excellent representation of the 
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structure factor near the convective instability. In this approximation: 



G NE (r||,L/2,L/2) = 
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where fy = r/L. The integral (52) can be evaluated exactly, when i? < i? c , 
but the result is a rather long expression. Instead we prefer to plot in Fig. 6 
the resulting Gne(?"||,z — z' — L/2)/Gne(0, £ — z' — L/2) as a function 
of r|| for two values of the Rayleigh number. The solid curve represents the 
density autocorrelation function for e = —0.01, which is close to the convective 
instability; the dotted curve represents the same function, but for e = —0.7. 
For R well below R c , the correlation function is overdamped but, as the value 
of R approaches R c , it starts to oscillate and the number of oscillations before 
the correlation function decays to zero increases the closer one approaches 
R c . To see this mathematically, we have made an asymptotic expansion of 
Gne^m, z — z' — L/2) for large r» and, when R is close to R c , we obtain: 
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The approximation (53) for e = —0.01 is shown in Fig. 6 as a dashed curve, 
to be compared with the solid curve representing the autocorrelation function 
as given by Eq. (52) for the same value of e. As can be inferred from Fig. 6, 
the approximation is excellent for small e and for large f\\. From Eq. (53), the 
real-space behavior of the autocorrelation function in the horizontal plane as 
a kind of damped oscillation is evident. 



8 Summary and concluding remarks 



Starting from the linearized fluctuating Boussinesq equations we have derived 
the structure factor of a horizontal fluid layer taking into account both gravity 
and finite-size effcts. The resulting expression reproduces the q~ 4 dependence 
of the nonequilibrium structure factor predicted theoretically [6-8] and con- 
firmed experimentally [10-13] for negative Rayleigh numbers, it accounts for 
the saturation of the nonequilibrium enhancement of the intensity of the fluc- 
tuations at small wave numbers observed by Giglio and coworkers [22,23,52,53] 
and it is consistent with the experimental observation by Ahlers and cowork- 
ers [28,51] of the structure factor close to the convective instability. We have 
thus provided a unified approach for describing nonequilibrium fluctuations 
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for both negative and positive values of the Rayleigh number, provided that 
R < R c . Our solution also enables us to evaluate the real-space behavior of the 
correlation function of the temperature fluctuations which agrees with the re- 
sults obtained by other investigators for the spatial direction co-incident with 
the direction of gravity [14,17,19,60,61]. 

There are still a number of issues that deserve further attention. 

First of all, for the sake of mathematical simplicity, we have adopted stress- 
free boundary conditions in the present paper. For a quantitative comparison 
with experimental data at very small wave numbers and close to the connec- 
tive instability modifications due to no-slip boundary conditions should be 
considered [25,27,29,31,36]. 

Second, we want to emphasize that neither our solution nor that of Swift 
and Hohenberg (SH) close to the convective instability yield a Lorentzian 
wave number dependence on q as is often assumed in the literature since a 
Lorentzian does not yield the proper asymptotic behavior for either small 
or large q. Moreover, even when the SH equation is not approximated by a 
Lorentzian, it does not yield the expected asymptotic q behavior either. We 
plan to study this issue in more detail in the future. 

Finally, while the average local values (9) and (w) are zero, the average cross 
correlation (6 ■ w) outside equilibrium is not generally zero. As a result, the 
coupling between temperature and transverse- velocity fluctuations may cause 
a small convective heat transfer, even below the instability [29,31]. A satisfac- 
tory solution of this effect has not yet been obtained [30-33]. 
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